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Abstract 

The Fermat principle is used to define trajectories in nonhomogenous optical 
media. The Poincare model of the Lobacevskii geometry is derived. The index of 
refraction is determined for the light confined in the circular trajectory in the optical 
medium. 
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1 Introduction 

Trajectories of elementary particles are the basic ingredients of physics 
of elementary particles and cosmical rays. No elementary particle can 
exist without its trajectory. While in particle physics the trajectories of 
particles are determined by their parameters as mass, charge, spin, velocity 
and by the influence of the magnetic and electric fields on its motion, the 
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trajectories of light in geometrical optics are determined by the index of 
refraction of the optical media. 

Geometrical optics considers the propagation of waves of light, as the 
propagation of rays, completely divorced from their wave properties. In 
other words, geometrical optics corresponds to the limiting case of small 
wavelength, A — > 0. 

The fundamental equations of geometrical optics for the direction of the 
rays are derived for instance by Landau et al. (1988, 1982) for any quantity 
/ describing the field of the wave (any component of E or H). For a plane 
monochromatic wave, we have for / the form 

/ = ae t{k - r -" t+a) = ae l( - krXr+a \ (1) 

(we omit the Re as the real part of all expressions). We write the expression 
for the field in the form 

/ = ae't (2) 

In case the wave is not plane, but geometrical optics is applicable, the 
amplitude a is, generally speaking, a function of the coordinates and time, 
and the phase ip which is called the eikonal, does not have a simple form, 
as in (1). It is essential, however, that ip be a large quantity. This is clear 
immediately from the fact that eometrical optics corresponds to the limit 
A^O. 

Over small space regions and time intervals the eikonal ip can be 
expanded in series. The time derivative of ip gives the frequency of the 
wave: 

and the space derivatives give the wave vector 

grad ip = k (4) 

and consequently the direction of the ray through any point in space. For 
a steady monochromatic wave, the frequency is a constant and the time 
dependence of the eikonal is given by a term —cot. We then introduce a 
new function ip\ (also called the eikonal), such that 

if> = -ut + (u/c)if>i(x,y,z). (5) 
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Then fa is a function of the coordinates only, and its gradient is 



grad fa = n, (6) 

where n is a vector such that 

k = wn. (7) 

The magnitude of n is equal to the refractive index n of the medium. 
Hence the equation for the eikonal in ray propagation in a medium with 
refractive index n(x, y, z) (a given function of the coordinates) is 

The equation of ray propagation in a steady state can also be derived 
from Fermat's principle, according to which the integral 

f-B r B 

fa = / n.dl = / ndl (9) 

J A J A 

along the path of the ray between two given points A and B has a value 
less than for any other path between A and B. 

Equating to zero the variation of this integral, we have 

5 fa = [ B (5n.dl + n5dl) = 0. (10) 

J A. 

Let (!)r be a displacement of the ray path under the variation. Then 
Sn = ^r.grad n, 5dl = l.dSr, where 1 is a unit vector tangential to the 
ray. Substituting in 5 fa and integrating by parts in the second term (using 
the fact that = at A and B), we have 

rB rB 

Sfa = / ^r.grad n dl + / rii.ddr = (11) 

J A. J A. 

j B A ( grad „ _ M . 6rdl (12) 



Hence 



grad n.dl = (13) 



Expanding the derivative and putting dn/dl = l.grad n, we obtain 



3 



d\. 1 

— = — [grad n — l(l.grad n)\. (14) 
dl n 

This is the equation giving the form of the rays. 

We know from differential geometry that the derivative dl/dl along the 
ray is equal to N/i?, where N is the unit vector along the principal normal 
and R the radius of curvature. Taking the scalar product of both sides of 
(14) with N, and using the fact that N and 1 are Perpendicular, we have 

1 N gradn 

R n 

The rays are therefore bent in the direction of increasing refractive index. 
The velocity of propagation of rays in geometrical optics is in the 
direction of 1 and is given by the derivative 

u = du/dk. (16) 

This is also called the group velocity, the ratio uj/k being called the phase 
velocity. However, the latter is not the velocity of physical propagation of 
any quantity. 



2 Light trajectory from the Fermat principle 

The point A and B can be joined by the infinite number of lines passing 
from A to B. There is the shortest distance between A and B forming 
the segment of the straight line passing through A and B (Hilbert, 1902). 
The shortest distance between point A and B can be physically realized 
by the flexible but non-elastic fibre. The segment AB of a straight line 
can be prolongated in the direction AB, or BA in order to generate the 
straight line. The shortest trajectory y = y{x) between two point A(xi, y\) 
and B(x2,y2) in the Euclidean plane x-y is y = ax + 6, a, b being some 
constants which can be find as the solution of the Bernoulli izoperimetric 
problem of the variational calculus with functional 

F=r[l + y y'M X ; y'= d J (17) 
Jxi dx 

after insertion of it into the Euler-Lagrange equation 

dy dx \dif) ^ 
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and after its solution (Lavrentjev et al., 1950). 

The geodetic line y = y(x),z = z{x) from point A{x\,yi,z-\) to point 
B(x2, y% Z2) on the surface (p(x, y, z) = is the solution of the izoperimetric 
problem with the functional 

F = £{ii + y ° + z T 2 -Mx,y,z)}dx; v' = %*' = % ("J 

where A is the Lagrange multiplicator. 

The Fermat optical theorem states that the trajectory of light from 
point A to B in the optical medium is of the shortest optical length. At 
the same time the trajectory of the optical ray from point A to point B 
with reflection on the mirror in point C is of the shortest optical length in 
optical medium. 

The trajectory of light passing from A(x\, y\) to B(x2, 2/2) is determined 
by the Fermat principle, which states that the time from A(x\,yi) to 
B(x2, yi) is the result of the minimization of the functional 

T(y,y') = p ±- = f'^f dx , (20) 

where v(y) is the velocity of light. 

The functional T(y, y') is the solution of the Euler-Lagrange equations 

(18). 

If v = Ay, the solution of eq. (18) is in the form of circles forming the 
Poincare model of the Lobacevskii geometry: 

(x-C) 2 + y 2 = r 2 . (21) 

The last equation can be also derived for velocity v = Ax. Then the 
functional (20) is of the form 



tv >\ f X2 ds f X2 Vi + y' 2 , / 00 x 
T(y,y)= ~r^ = — i-^dx, 22 

J xi V(X) J xi viX) 



v(x) Jx\ v(x) 

We get from the Euler equation with y' = dy/dx 



T y - ^-T y , = (23) 



we get T v ' = const = C\, or 
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x^/l + y 

We get with the substitution t = arctany', 



y' 2 

It follows from ?/ = tan/; that % = C2 sintdt. Or 



x = — — H . = — sin t = C2 sin t. (25) 



2/= -C 2 cos£ + C 3 . (26) 

Or, 

x = C2sin£, y — C3 = — C2 cost. (27) 
If we eliminate the variable t in the last equation, we get 

x 2 + (y - C 3 ) 2 = Cl = r\ (28) 
which is a circle with the center on the y-axis. 

3 Trajectory of light in the stratified medium 

Let us consider the stratified 2D medium in the plane x — y : where every 
layer is parallel with the x-axis. Then, dx — > (Ax)i : dy — > H/n, in (20) 
where i = 1, 2, 3, and H is the height of the stratified medium. In this 
case the velocity is Vi = v(yo + iH/n). Then, instead of functional T(y, y') 
we write the following sum: 



T _ g y/((Ax)i)* + (Ay) 2 _ g ^(a l+l - atf + (H/n) 2 ^ 

i=0 v i i=0 v i 

The last sum has the stationary value if and only if dT/dcii = 0. Or, 

dT/da, = 



(Oi+i - Oj) 



(cii - cii-i) 



Vi\f(a i+ i - di) 2 + (H/n) 2 Vi-iJ(a,i - a f _i) 2 + (H/n) 2 



COS Ifj COS (£j-i _ Q 



(30) 



Or, 



cos (pi cosy?j_i I 

= const = — . (31) 



Vi Vi_i ' k 

Or, in the continual limit, 

cos (f 1 
v(y) k 

However, for y = y(x) it is tan<£> = y' and cose/? = -j= n 
Then we have 



X - X °-J7W=^- (33) 



We get for v = Ay: 



'j/o yjk 2 — A' 2 y 2 



with the solution 



(32) 



( X -Cf + i={^ = R\ (35) 

which is the circle with the centre on the x-axis. 
For v = A/y, we get instead of equation (34): 

r Ady 

X = J VkV-^ + (36) 

Let us remark that the physical meaning of the relation (36) with 
v = A/y is the problem of the catenary in the homogeneous gravitational 
field. 



4 The Poincare optical model of the Lobacevskii 
geometry 

The Lobacevskii geometry is the integral part of the general geometry 
called non-euclidean geometry. The name non-Euclidean was used by 
Gauss to describe a system of geometry which differs from Euclid's in 
its properties of parallelism. Such a system was developed independently 
by Bolyai in Hungary and Lobacevskii in Russia, about 120 years ago. 
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Another system, differing more radically from Euclid's, was suggested 
later by Riemann in Germany and Schlafli in Switzerland. The subject 
was unified in 1871 by Klein, who gave the names parabolic, hyperbolic, 
and elliptic to the respective systems of Euclid, Bolyai-Lobacevskii, and 
Riemann- Shlafli (Coxeter, 1998). 

The substantial mathematical object in the Lobacevskii geometry is the 
angle of parallelism defined by Lobacevskii as follows. Given a point P 
and a line q. The Intersection of the perpendicular through P let be Q 
and PQ = x. The intersection of line p passing through P, with q, let be 
R and QR = k. Then, the angle RPQ for perpendicular distance x 

U(x) = 2tan- 1 e- x/k . (37) 

is known as the Lobacevskii formula for the angle of parallelism (Coxeter, 
1998; Lobacevskii, 1914). 

The Poincare model of the Lobacevskii geometry is the physical model 
of the optical trajectories in a medium with the velocity of light v = Ay. 

According to Hilbert (Hilbert, 1903), it is not possible to realize the 
Lobacevskii geometry globally on surface with the constant negative curva- 
ture. The Beltrami realization of the Lobacevskii geometry is only partial. 
The Lobacevskii geometry is the partial geometry on the pseudosphere 
with the parametric equations 

x = a sin w cos y = a sin usinv, z = a ^lntan ^ + cos uj . (38) 

(Kagan, 1947, ibid. 1948; Efimov, 2004; Norden, 1956; Klein, 2004; 
Manning, 1963). 

The pseudosphere is the surface generated by the rotation of the Leibniz 
tractrix with equation 

x = asinu, y = 0, z = a ^lntan ^ + cos iij . (39) 

The pseudosphere is in the half geodetic coordinates given by the 
squared element 

ds 2 = du 2 + cosh 2 -dv 2 . (40) 

a 

The pseudosphere is in the izothermical coordinates given by the 
Poincare squared element 
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2 a(dx 2 + dy 2 ) 



(41) 



The Leibniz solution of the tractrix problem is as follows: 



y = a 



ln(a + Va 2 — x 2 ) 



— \/ a 2 — x 2 . 



(42) 



x 



The trajectory of light in the Poincare model is a trajectory passing from 
A{x\,yi) to B(x2,y2) and determined by the minimal time from A(x\,yi) 
to B{x<i,y<i). It is the result of the minimum of the functional (20, 22). 
The Poincare circles in his model are analogue of the for straight lines in 
the Euclidean geometry. The following theorems are valid in the Poincare 
model of the Lobacevskii geometry: 

Theorem 1:. Only one half-circle passes through two points A, B in 
the Poincare plane. 

Theorem 2: The curvilinear segment AB in the Poincare plane is of 
the shortest length. 

Theorem 3: The parallels are two half-circles with the intersections on 
the x-axis. 

Theorem 4: If point A ^ q then there are q\ \\ q, q 2 || q passing through 
A, with qi ^ q 2 . 

Theorem 5: If point A q, q\ \\ q, q 2 \\ q, then qi, q 2 divide the Poincare 
plane in four different sectors I, II, II, IV. 

Let us remark that he optical distance between point A and B is not 
equivalent to the mechanical distance realized by the nonelastic flexible 
fibre as the shortest distance between point A and B. The Poincare model 
of geometry where the light velocity is v = Ay is the interaction model of 
light with the optical medium. 

It is elementary to see that if we define the Poincare problem on a 
sphere, then we get so called spherical Poincare model of the Lobaevskii 
geometry. 

Let us use the index of refraction n(r) in the Euclidean plane and polar 
coordinates r, if to derive the Poincare model. 
The the explicit form of the Fermat principle 




(43) 



is (Marklund et al, 2002) 
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5 J n(r) 



\ 



MS) =o 



(44) 



The last equation is equivalent to the Euler-Lagrange variational equa- 
tion for the functional F(ip, tp') 



K - ^ - o. 



(45) 



Or, 



d 
dr 



n(r) 



r 2 dipdr 



= 0; 



/ 1 + ^ 2 ©' 

It is evident that the elimination of d(p/dr is as follows: 

C 



(46) 



dr 



= 0. 



(47) 



The circular trajectory is defined by equation dr/dtp = from which 
follows the index of refraction as 



n(r) = 



const 



(48) 



The experiments has been performed in Bose-Einstein condensate with 
the result that the optical light pulses can travel with extremely small 
group velocity about 17 meters per second (Hau et all. 1999). 



5 Discussion 



We explained the Fermat principle in the geometrical optics and in the 
variational calculus. We determined the optical trajectories in the optical 
media. We also derived the Poincare model of the Lobacevskii geometry. 
The Poincare model of the Lobacevskii geometry is the optical model based 
on the index of refraction, which is the consequence of the interaction 
of light with medium. Also, space, when considered independently of 
measuring instruments, has neither metric nor projective properties; it has 
only topological properties. It is amorphous (Poincare, 1963). That which 
cannot be measured cannot be an object of science. We applied the Fermat 
principle to the determination of the index of refraction of the medium to 
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get the circular optical trajectories. This effect is the confinement of light 
by optical medium. 

There is no doubt that the monochromatic optical beam is composed 
from photons of energy E = hu. While the rest mass of photon is zero, 
the relativistic mass follows from the Einstein relation E = mc 2 . After 
identifying the relativity energy and quantum energy of photon we have 

hco 

m = 49 

& 

The centrifugal force acting on photon moving with velocity v in optical 
medium along the circle with radius R is for the photon mass as follows: 



■^centrifugal ^ ~E> ' 



flUJ V' 

The centrifugal force is the origin of the unstability of the photon 
trajectory in the optical medium. The further origin of the unstability 
of the trajectory are the thermal fluctuations of the index of refraction. 
We know, that the Kapitza effect is based on the thermal fluctuations 
of the index of refraction (Landau, et al., 1982). So, the experimental 
investigation of the confinement of photon in the optical medium is 
meaningful at T « 0. There is no doubt that the investigation of the 
photon trajectories is the crucial problem of the optical physics and it is 
interesting for all optical laboratories over the world. 
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